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The thermal conductivity of the (2+l)-dimensional NJL model in the presence of a constant 
magnetic field is calculated in the mean-field approximation and its different asymptotic regimes 
are analyzed. Taking into account the dynamical generation of a fermion mass due to the mag- 
netic catalysis phenomenon, it is shown that for certain relations among the theory's parameters 
(particle width, temperature and magnetic field) , the profile of the thermal conductivity versus the 
' applied field exhibits kink- and plateau-like behaviors. We point out possible applications to planar 

O ' condensed matter. 

■ 

1 ■ PACS numbers: ll.lOKk, ll.SOQc, ll.lOWx, 12.20.Ds 
■ 
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It is now well established, from the study of many relativistic theories of massless fermions in the presence of an 
external magnetic field, that a magnetic field can be a strong catalyst for chiral symmetry breaking with the consequent 
generation of a fermion dynamical mass even at the weakest attractive interaction among fermions . This magnetic 
^ , catalysis (MC) of chiral symmetry breaking has proven to be universal, its main features being independent of the 
Q model under consideration (see Refs. 0— § for various aspects of this phenomenon). The universality character 
, of the magnetic catalysis has motivated many recent works JT3| aimed to apply it to diverse areas of quantum 
physics. The essence of the MC effect lies in the dimensional reduction of the fermion pairing dynamics due to the 
' confinement of these particles to their lowest Landau level (LLL), when the pairing energy is much less than the 
J> Landau gap ^/B {B is the magnitude of the magnetic field induction). Under these circumstances, any attraction 
\ between fermions, whenever small it might be, is strengthened by the effective dimensional reduction in the presence 
' ' ' of the magnetic field, and therefore, a condensate of fermion- antifermion is formed with the subsequent generation of 
ijj a fermion mass. The lowest LL plays in this case a role similar to that of the Fermi surface in BCS superconductivity 

■ It is known that several quasiplanar systems have low-energy excitation spectrum of quasiparticles (QP), charac- 
^3 ■ terized by a linear dispersion, around the Fermi surface consisting of isolated points. The dynamics of these QP can 

■ be described by a "relativistic" quantum field theory of massless fermions. When such Dirac-like QP are electrically 
\ charged, they can couple to an externally applied magnetic field which can catalyze the condensation of QP-antiQP 

H • pairs. Then, one would expect the realization of MC in such a kind of condensed matter systems. As a matter 
\^ ' of fact, the MC was suggested as the possible explanation jl^] for the profile of the thermal conductivity in 
an applied magnetic field observed in recent experiments in planar high-Tc cuprates [0. The MC has also 

been proposed jlj] as the source of the semimetal-insulator phase transition observed in the so-called highly oriented 
pyrolitic graphites (HOPG) in the presence of a magnetic field. 

Given that the heat transport is a convenient probe to understand many basic properties of quasi-planar condensed 
matter systems, as gap structure, QP density, scattering rate, etc, the study of heat transport in a quasiplanar system 
subjected to MC may turn out physically revealing. It is the goal of the present paper to study how the MC affects 
?-j \ the thermal transfer properties of a (2-|-l)-dimensional fermion system under an applied constant magnetic field. In 
particular, we calculate the thermal conductivity of a system of QP described by a (2-|-l)-dimensional Nambu-Jona- 
Lasinio (NJL) model which exhibits a mass (gap) generation for fermions in the presence of a magnetic field, and 
discuss possible applications to planar condensed matter systems. To obtain the results here reported, we used the 
same approximation of a constant magnetic field that was already explored in Refs. However, our calculations 

deviate considerably from what was done in these papers. Not only we take into account the contribution of all 
Landau levels, but the definition we use for the heat current itself is different. In our formulation, when the gap 
induced by the magnetic field is opened, the thermal conductivity exhibits a new term proportional to [a is the 
gap). Near the phase transition point, the gap behaves like a ^ y/ B ^ Be in the mean-field approximation. Hence, 
the term proportional to yields a positive contribution in the slope of the thermal conductivity, leading to a jump 
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in the slope of k{B) at B = Be (kink-like behavior). Notice that the magnetic catalysis is mainly responsible for the 
kink effect, being the critical behavior of the gap an essential factor for this result. 

We underline that to obtain such a kink-like behavior, it is crucial to go beyond the LLL contribution in the 
calculation of the thermal conductivity, since, as we show, the heat transfer takes place due to transitions between 
neighboring LLs. Physically, this is easy to understand keeping in mind that the spatial momentum of the QP lies 
in the system plane. In the presence of a magnetic field perpendicularly applied to the two-dimensional sample, the 
QP spatial momentum is purely transverse and hence quantized into discrete Landau levels. Therefore, the transfer 
of kinetic energy can only occur by means of transitions between Landau levels. 

Even though our results were obtained by using a particular model, we point out that the main outcome of the 
present work is of a more general and theoretical character, as we show that the MC phenomenon can be responsible for 
a kink-like effect in the thermal conductivity of a whole class of (2-|-l)-dimensional relativistic fermion systems. That 
is, we show that the kink effect is essentially model independent, since it is determined by the critical behavior of the 
dynamically generated mass near the phase transition point. This fact makes the basic outcome of our investigation 
relevant beyond the particular model under consideration, linking it to the universality class of theories with such a 
critical behavior. In connection with this we conjecture that, since the HOPG may be described by a model that 
belongs to the same universality class as that of the model used here, the HOPG thermal conductivity in the presence 
of a magnetic field should display similar kink-like behavior. 

The plan of the paper is as follows. In Sec. || we derive the expression for the thermal conductivity in the (2+1)- 
dimensional NJL model in the presence of a constant magnetic field and analytically study its different asymptotics. 



underlying the possible application of each result. In Sec. Ill we obtain the thermal conductivity vs magnetic field 
profile using numerical calculations. In the reported graph, the change of slope in the thermal conductivity profile is 
shown to occur at the critical magnetic field where the fermion dynamical mass is generated at the given temperature. 
The conclusions and discussion of potential applications of our results are presented in Sec. IV. In Appendix ^ we 
derive the critical curve in the B — T plane which separates the symmetric and the symmetry-broken phases in the 
(2-t-l)-dimensional NJL model, and study the scaling of the mass near the critical curve. A derivation of the Kubo 
formula in the framework of Matsubara formalism is given in Appendix B. 

II. THERMAL CONDUCTIVITY IN THE (2+l)-D NJL MODEL IN THE PRESENCE OF A CONSTANT 

MAGNETIC FIELD 

A. General Results 

We start from the (2+l)-dimcnsional NJL Lagrangian density in an external magnetic field 

/: = i[^„z7^z?^^.] + ^(^,V.)'. (1) 

where = — ieA'^^* is the covariant derivative and the vector potential for the external magnetic field is taken 
in the symmetric gauge 

^r- (o.-f^2,fxi). (2) 

We assume that the fermions carry an additional flavor index i — I, . . . , N (N = 2 for realistic d-wave superconduc- 
tors). The Dirac 7-matrices are taken in the reducible four-component representation. 

In the absence of the bare mass term rmpip, the Lagrangian density (|^) is invariant under discrete chiral symmetry 

V' 75^^, -0 -i'lb, (3) 

which forbids the fermion mass generation in perturbation theory. The appearance of the mass (energy gap) is due to 
the spontaneous breaking of the above discrete symmetry that leads to a neutral condensate of fermion-antifcrmion 
pairs. In condensed matter physics this could correspond to the condensation of excitons (electron- hole bound states). 
Introducing the composite field a = —g{'ijji^i)/N, the Lagrangian (|l|) can be written in the form 

£ = i [i,,,ij^^D^i,,] - a^.^H - (4) 

One readily verifies the equivalence of the Lagrangians (|^) and (^ by making use of the Euler-Lagrange equations (or 
performing the integration over the field a in the functional integral). The field a has no dynamics at the tree level. 
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however, it acquires a kinetic term due to fermion loops. The vacuum expectation value of a gives a dynamical mass 
(gap) to fermions. The effective action for the composite field a can be obtained by integrating over fermions in the 
path integral (see Appendix^. 

It is well known that in the absence of a magnetic field and at zero temperature the mass generation occurs only if 
the coupling constant exceeds some critical value pQl . This can be seen from the stationary equation for the effective 
potential corresponding to the Lagrangian density (|l|) , which in leading order of 1 /A^-expansion is given by 



dVjcj) 



A A 

- 

g 



(5) 



(it follows from Eq. (Q in the limit B = T = 0). In (|) we introduced the dimensionless coupling constant g — Ag/n 
and A is the ultraviolet cutoff parameter. From Eq. (H) it is easy to see that there exists a critical value for the 
coupling constant — such that, if g < g^, Eq. (0) has only the trivial solution (ct = 0), while for the strong 
coupling limit g > g^, a nontrivial solution (ct = A{'g — y/n)/ 'gy/ir) is reached which leads to the generation of a 
fermion mass. Hence, the critical coupling separates two phases: the weakly coupled massless phase atg<g^ and 
the strongly coupled massive one (g > ffc)- ^.n applied magnetic field changes the situation dramatically, so that the 
mass generation now takes place at all g > |jl|-^|, hence the name magnetic catalysis. At finite T and B ^ there 
is a critical curve in the B — T plane separating the symmetric and the symmetry broken phases (the derivation and 
analysis of the critical curve in the B — T plane are given in Appendix |A|). 

To derive an expression for the static thermal conductivity in an isotropic system we follow the familiar linear 
response method and apply Kubo's formula 

oo 

K=-^^Im j dtt J (fxicfx2{ui{xi,0)ui{x2,t)), (6) 



where V is the volume of the system, T — is the temperature, and Ui{x,t) is the heat-current density operator. 
The brackets denote averaging in the canonical ensemble with the density matrix p = e~^^ /Z, Z = Tie~^^ . 

Physically the thermal conductivity k appears as a coefficient in the equation relating the heat current to the 
temperature gradient 

u = -kVT (7) 

under the condition of absence of particle flow. If we neglect the chemical potential the heat density coincides with 
the energy density, hence the quantity that satisfies the continuity equation 

e{x) + V-u(a;) = (8) 

can be interpreted as the heat current density. Equation (^) defines u to within a divergenceless vector, which is 
sufRcient for calculating the conductivity. The vector u is obtained automatically from the Lagrangian density (|l|) as 

dC ■ - dC i - 

^ + = o (^7^5o^ - doM) ■ (9) 



9(a*V) d{d'^) 2 
When using equations of motion the last expression can be represented in the form 

u,^'- {^j°D,i, - A?7V) ■ (10) 

At this point it is useful to underline that our definition of the heat current does not coincide with the one used in 
Rcfs. Pi{x) = ^pj^diijj — ditpj^ip. Operator Pi{x) cannot be obtained from Equation (||) by using the equations 

of motion, unless the external field is zero, so it does not lead to the correct thermal conductivity in the presence of 
a magnetic field. Note also, that our quantity Ui is explicitly gauge invariant in contrast to Pi. 

The correlator of the heat currents (or polarization function) in (^ is evaluated in the following way (see Appendix 
B for details on the formalism). First, it is computed in the Matsubara finite temperature formalism replacing the 
time t by the imaginary time t [t — —it): 

P 

= ly" dTe*""-(r,;7,(T){/,(0)), c/,(t) = j d^xu,{x,T), = (ii) 
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The thermal conductivity is then given by the discontinuity of the retarded function n^(ri), which is obtained by 
analytic continuation from imaginary discrete frequencies 11^(17) = n(if2m — > f2 + ie): 



(12) 



Neglecting vertex corrections]^, the calculation of the thermal conductivity reduces to the evaluation of the bubble 
diagram 



-tr 



(13) 



where S{uj, k) is the Fourier transform of the translation invariant part S{x — y) of the fermion propagator in an 
external magnetic field: 



S{x - y) = exp ie / S{x - y). 



(14) 



Note, that the translation non-invariant phase of the fermion Green's function cancels in the computation of H. 
Defining a spectral representation for 5(iw„, fc), we can write 



duj A(uj, k) 



lUJn — ^ 



(15) 



The spectral representation allows one to make analytic continuation and find the retarded, S^, and advanced, S^^ 
Green functions according to the rule S^'{ijJ + ie, k) = S{iuJn = oj + ie, k) and S^{ijJ — ie, k) — S{iuJn = lu ~ ie,k). The 
spectral function A{lo, k) is given by 



A{uj,k) 



2TTi 



S'^{uj-ie,k)-S^{uj + ie,k) 



-ImS^iuj + ie,k). 



(16) 



Plugging the spectral representation ( [l5| ) into Eq. (|T^) the sum over Matsubara frequencies can be performed 0. After 
this has been done, we can continue the external frequencies to the real axis to get 11^(17). Finally, we arrive at the 
following expression for the thermal conductivity 
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327rT2 J cosh^ 



2T 



d^fctr 



YA{uj,k)j'A{Lo,k) 



(17) 



In order to compute k we have to specify now the fermion propagator. Since the fermion mass is generated in the 
(2-|-l)-dimensional NJL model in a magnetic field already at weak coupling, we take the standard expression for the 
massive fermion propagator in a magnetic field, decomposed over the Landau level poles P,p3[ 



^-^ uj^ — (y^ — 2eBn 

n— 



(18) 



where 



^It has been argued that for small impurity densities the thermal conductivity, unlike the electric conductivity, is unaffected 
by vertex corrections |Q . 

^There is a subtle point here since the sum over frequencies appears to be divergent. However, as was shown in this 
divergence results from an improper treatment of time derivatives inside the time-ordered product of currents in (|ll|)T This 
divergence disappears when the problem is treated more carefully. The prescription is simply to ignore it. 
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2fc2 



eB 



eB 



(19) 



with P± = (1 ± 17^^7^)72 being projectors and i„,L^ Laguerre's polynomials (iLi = 0). Here a is the ferniion 
dynamical mass obtained from the finite temperature gap equation in a constant magnetic field (see Appendix ^) . 
The spectral function according to (n6h is found to be 



A(w,fc) 

where M„ = Vct^ + 2eBn and 



_ii r ^ (-1)" 

n— 



(7°M„ + a)/i(fc) + /2(fc) , (7"M„-a)/i(fc)-/2(^) 



/i(A:) = 2 



2fc2 



(a; - M„)2 + n 



2fc2 



P_L„ — - P+i„-i — 



eB 



eB 



(w + M„)2 + r2 



/2(fc) = 4A7i^_J^ 



(20) 



(21) 



Here, we introduced the width F of the quasiparticles, which is due to interaction processes, in particular, scattering 
on impurities, having replaced e in ( p^ ) by a finite T. In general, the scattering rate F, which is defined through 
the fermion self-energy, r{uj) — — Imi;^(a;), is a frequency-dependent quantity (as well as temperature and field 
dependent). It must be determined, like the dynamical mass, self-consistently from the Schwinger-Dyson equations. 
At low temperatures we are interested in its value at u; = 0, so we will consider it as a phenomenological parameter. 
In the absence of a magnetic field {B = 0) it can be shown that the spectral function (20) reduces to 



A{uj,k) 



T 
2^ 



"f"E -k^ + a "f"E + k^~ a 
{lu - E f + r2 ^ [uj + E f + r2 



E 



hence, the retarded fermion Green's function is 



S^{iu,k) = 



7o(w + iT) — k^ + a 
(w + iVy -k^-a^' 



Straightforward calculation of the trace in (|17|), with A{uj, k) from Eqs 



(22) 



(23) 



gives 



tr 



-f'A{uj,k)fA{uj,k)\ = 



2fc2 
eB 



n,m— 

^ 2^2' 



m+n+1 



(c^2 + M2 + r2)(a;2 +MI+ r2) - 4C^2^2 
[(^2 + M2 + r2)2 _ 4lj2M2][(w2 + ^1 + T^f - AlO^ M^] 



Lm-1 



eB 



+ Ln-l 



2fc2 
'eB 



2fc2 



Performing now the integration over momenta in Eq. (|r^) produces Kronnecker's delta symbols Sn,m-i 
to the orthogonality of Laguerre's polynomials, thus we get 



eBT^N 



E 



(c^2 + Ml + r2)(c^2 ^ Ml^^ + r2) - 4C^2^2 



cosh 



2T 



[(^2 + Ml + F2)2 _ 4lo2M2][(c^2 + ^,^2^^ + r2)2 - 4^2^^^ 



(24) 
%-i due 

(25) 



n+lJ 



Note that the factor eB in front of the right hand side of ( psj) originated from integrating over transverse momenta and 
gives the degeneracy of Landau levels (more exactly, the degeneracy is NeB/2'iT for the lowest LL, n = 0, and NeB/ir 
for levels with n > 1). We stress the important result that because of the appearance of the mentioned Kronnecker 
deltas only transitions between neighboring Landau levels contribute into the heat transfer. Had we restrict ourselves 
to the LLL as in we would have gotten zero result for k. 

Further summation over n in Eq. ( p5| ) can be performed expanding the integrand in terms of partial fractions. The 
resulting sums are expressed through digamma functions by means of the formula 



E 

n=0 



A 
n + a 



B 



C 

n + c 



D 



= - [A^j{a) + Bi>{b) + Ci>{c) + Dij{d)] , 



(26) 
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where for convergency A + B + C + D = 0. 

After some algebraic manipulations the final expression for k is written as follows 



2tt^T^J cosh^ ^ (eS)2 + (2c^r)2 



+^2^ r2)(eB)2 - 2uj^{oj^ -a^+ T^)eB 



(w2_cr2-r2)2 + 4t^2r2 



C^(cj2-a2^p2) 



2iujT 



2eB 



(27) 



This formula is the main result of our paper. Note that it is independent of the particular model used to describe 
the QP interactions unless we specify the dependence of the dynamical mass on T,T, eB for a concrete model. 

Another representation of Eq. ( |27| ) which is particularly convenient for studying the small width limit T <^ T, \/eB 

is 



Arr2 



1 



(eBf 



^ M 



cosh 



2T 



(eS)2 + {2ujTy 



+ eBuj{uj + cr) + eBuj{uj - a) 



(uj + (t)2 + r2 



{uj - ct)2 + r2 



n=l 



M2 



■ 2wM„ 2LuMn - 0-2 _ ^^2 



(w + M„)2 + r2 (co-M„)2 + r2 



(28) 



It is obtained from ( p7| ) if one takes the series representation of the ?/;-function in the integrand of (^) and writes the 
expression in the curved brackets in fractions of l/(r2 + x'^). 

We are now in the position to study different asymptotic regimes defined by different relations among the dimensional 
parameters a,T,T, B. 



B. Zero Magnetic Field 

First, we consider the limit of vanishing magnetic field (-6 = 0). For that, one can use the formula for the asymptotic 
of ■i/'-function at large values of the argument 



'ip{z) = log z ■ 



1 



1 



1 



2z 12z2 



120z4 



O 



1 



(29) 



to get 



A^ 



An^T^J cosh^ 



2T 



2wr 



vr a2 ^ p2 _ ^2 

arctan 

2 2wr 



(30) 



This expression describes the behavior of the thermal conductivity as a function of temperature for massive Dirac's 
particles and is relevant for the supercritical phase of the NJL model, where the mass is generated spontaneously even 
at zero magnetic field. 

The last expression can be evaluated analytically in two regimes, F <C T: 



A^ 



87rr2F 



duj- 



to (co^ - a^) 



cosh2 {lu/2T) 



-e T ^ T « a, 



(31) 



and F > T: 



Kq^N 

T ~ 3 



7^2 2^2p2(p2^5^2) 



(32) 



_CT2 + r2 15 (o-2 + r2)3 

Accordingly, in the weak coupling phase of the NJL model where the dynamical mass is not generated we obtain: 



T TT 



9C(3)T ln2F 
4 F ^ ~Yt 



F<r, 



(33) 
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and 



T 3 



^ 15 r2 



r > T. (34) 



Eqs. (p3|), ( p4| ) up to an overall factor coincide with the corresponding expressions obtained in Ref. [|4j for the 
vortex state of nodal quasiparticles in the d— wave superconducting phase of high-Tc cuprates. The overall factor 
there equals (v\, + v\)/vpV/\ where vp,v/^ are respectively the velocities perpendicular and tangential to the Fermi 
surface. They originate from the quasiparticle excitation spectrum in the vicinity of the gap nodes which takes the 
form of an anisotropic Dirac cone E{k) — \/vpkl + f^fc|. With the overall factor replacing (= 2 in real d— wave 
superconductor), the first term in Eq. (|3j) reproduces the universal (or residual) thermal conductivity at low T in 
the so-called "dirty" limit, T <^ T, since it is independent of the impurity density, thus it will not depend on the 
specific characteristics of the scattering processes in the sample. The residual conductivity was recently observed in 
experiments p5[ | confirming the existence of g apless quasiparticles in c?-wave cuprates at T < Tc- Note, however, that 
in contrast to what was claimed in Ref. pq| , the low-temperature thermal conductivity for massive quasiparticles 
(Eq.(^2|)) does not exhibit a universal behavior when T — > 0. This peculiarity of the low-temperature thermal 
conductivity can be used to find out experimentally the second gap in cuprates. 

Expressions ( |33| ) , (^4|) were recently used to propose a scenario for the arising of a plateaux at high magnetic fields 
in Krishana's experiment. In that scenario the width F of QPs becomes dependent on the field due to scattering 
on disordered vortices, thus F becomes Fq -I- Tb where the field induced width F^ is calculated to be ^ ^/B- All 
the information regarding the magnetic field is encoded now in the total width F, hence the magnetic field is not 
explicitly present. If we start with a weak magnetic field , when F^ ^ Fq <C T , the thermal conductivity follows 
first the expression (|3^) (weak field regime) decreasing with the field. At some point, when F becomes of order T a 
crossover takes place to the high field regime (^) with plateaux. Physically, such a scenario is applicable only if there 
is a small number of vortices with large distances between them and the magnetic field is basically confined in tubes. 
However, for the field range of interest, Hd <^ B <^ Hc2, where Hci,Hc2 are the lower and upper critical magnetic 
fields of the high Tc superconductor respectively, the vortices are dense enough to overlap strongly giving rise to an 
effective uniform magnetic field in the whole plane ||27|| , so the above scenario is not perhaps the most appropriate for 
this field range. 



C. Non-zero Magnetic Field 



We shall analyze now the thermal conductivity in the presence of a uniform magnetic field in the whole plane. The 
analysis in all cases will be made at fixed T and F, and we do not assume the last one to be dependent on the field. 

From the phase transition analysis of the (2-|-l)-dimensional NJL model (see Appendix]^ it follows that at finite 
temperature there exists a critical value of the magnetic field Bc{T), above which the magnetic catalysis phenomenon 
occurs generating a dynamical fermion mass even at weak coupling (in what follows we consider only the weak coupling 
case g < gc)- For magnetic fields less than the critical one {eBc{T) ^ IQT^) the dynamical mass is zero (cr = 0). 



1, Narrow Width 

To study the narrow width limit F — > 0, we replace the fractions F/(F2 + x^) in Eq.(|2^) by t:5{x). Then, after 
integrating over lj, what is equivalent to evaluating in the mass shell for the different LL, we obtain 

^ NT { {eBf 0-2 ^ (2eS)2n(cr2 + 2eBn) 1 

In (|35| ) we kept F^ in the denominators in order to be able to reproduce a smooth behavior of k{B) in the limit 
B ^ 0. The origin of the first term in (|35|), whose essential role in the kink-like behavior is discussed below, can 
be traced back to the leading contribution of the zeroth to first LL transitions. Note that it contributes only when 
the dynamical mass is present, i.e. if a ^ . That is, because the ratio F/((a; — (t)^ + F^) appearing in Eq.([27|) 
becomes TrS{uj — a) in this limit, it does not contribute to ( [3^ ) unless cr ^ 0, due to the presence of the factor wMn 
the integrand of Eq. ( p7| ) . This means that in the narrow width limit the magnetic catalysis is not only connected to 
the generation of the mass, but it is responsible also for the enhancement of the transitions between zeroth and first 
LLs. The fact that the thermal conductivity is proportional in this limit to the scattering rate (width) F means that 
it results from transitions of quasiparticles between cyclotron orbits mediated by scattering of QPs on impurities. 
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a. Weak Field Limit, V eB < 4T. 

Taking into account that at weak couphng no dynamical mass is generated for fields below the critical value 
{eB < eBc), we take cr = in the calculation that follows. Making use of the Euler-MacLaurin formula 

oo 

-F(0) + Fin) c / dxF{x) - -F'{0), (36) 



expression ( p5| ) with zero dynamical mass can be recast for subcritical fields V eB < \/eBc ^ 4r in the form 

C30 

2NT^ r dxx^ 1 



7rr J cosh a; t2 + (-£^^ 



(37) 



Eq.(37) shows monotonic decreasing of n with increasing magnetic field B (at i? = it reproduces the leading term in 
Eq.(33)). Note that the scale y/ATT marks the crossover point where the transition from superweak (V eB < \/4TT) 
to weak (ViTT < < 4r) fields takes place. 

b. Strong Field Limit, \J eB > AT . 

We shall consider now the strong field regime, \feB > 4T, where a nonzero dynamical fermion mass is generated in 
the weakly interacting system (we are interested mainly in the region of coupling constants g 5c where the scaling 
a ~ \feB is achieved) . 

Let us start, however, analyzing the case of free massless fermions (ti = 0). In this case, one can use Eq. (|35|), after 
evaluating it in cr = 0, to describe the very large field (%/ eB ^ 4T) behavior of k in the narrow width case, what 
yields an exponential fall of the conductivity 

^NVeB 

K ~ —^e T . (38 

Coming back to the interacting case and after dropping the term depending on T in the denominators of Eq.(p5|), we 
obtain 

4^r2 I cosh^ ^ cosh^ J ^ ' 

In the limit of large fields (V eB ^ 4r), the first term in ( p9| ) is the leading one. Such a term would produce a sharp 
plateau, were the fermion mass a constant. Since in the model under consideration the mass is dynamical and it 
behaves a.s a ^ \feB at \feB ^ y/eBc — 4T, the first term gives rise to an exponential decrease as in the case of free 
massless fermions (^8|) for asymptotically large fields. 

c. Near the phase transition point, eB > eB^- 

In this case, the thermal conductivity can still be approximated by Eq.(p9|). As shown in Appendix A, near the mean 
field phase transition point a « ^^/eB — eBc, so if the field lies in the interval eBc < < 2eBc, the dynamical mass 
a < 2T and hence the cosh^ ^ appearing in the first term of Eq.(|39|) is of order one. In this field region the first term 



of Eq.(39) gives a positive contribution to the derivative of k close to the critical point. That positive contribution 



leads to a jump in the slope of k at eB ~ eBc therefore showing a kink-like behavior for k in the narrow-width case. 



2. Finite Width 



Let us consider now the case where the width F is small but finite. We are particularly interested in the behavior 
of K near the phase transition point, where eB > eBc, and therefore V eB > AT, F, cr. From Eq. (p5|) one can see 
that for these fields the contribution of transitions between Landau levels with n > 1 in the integrand behaves as 
l/(ei?)^,while the transitions between zeroth and first LL decrease as ^ 1/eB. However, the LL degeneracy is also 
proportional to eB, what implies that the transitions between the zeroth and first LL are not suppressed despite the 
fact that the gap between levels grows with the field. The leading in 1/eB behavior is easy to obtain from Eq. (p8|) . 
It is given by the expression 

C30 

NT^ r dujLu^ r/ 2lu{uj + (7)\ 1 / 2uj{uj-a)\ 1 

" ~ 4^ J cosh^ I V ^ ^ J'C^T^W+^^V^ ^ ) - + r2 
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As discussed in the previous subsection, near the transition the dynamical mass behaves as a 
can expand Eq.(^o|) around cr = to obtain 



hy/eB — eBc, so we 



dujLu 
cosh^ ^ 



1 + 2Lj^/eB 



3io^ - 
tj2 + r2 



2u/ 

eB 



U) 



3r^ 



r2 



(41) 



Clearly, for T < ^/3{2T), the term proportional to gives a positive contribution to the derivative of k with respect 
to B near the transition point, so there is a jump in the slope of k at the critical point: a kink-like effect. Notice 
that if cr were zero or constant, the derivative of the thermal conductivity would satisfy ~ (eB)'^ ,with C positive, 

so no kink-like effect would be present. On the other hand, since a is dynamical, near the transition point ~ 
a — P/eB + 0{{eB)^'^), with a and /? positive, hence one can see that the dynamical mass not only allows for a jump 
in the slope, but it flattens the profile after the critical field, allowing at least for certain region of the parameter 
space a behavior of almost zero slope: a plateau-like profile. This kink and plateau -like behavior is corroborated by 
numerical calculations in the next Section. 

We highlight that the dynamical mass is needed to obtain the kink-like effect in both narrow and finite width cases. 
Hence, in our model the kink of the thermal conductivity is directly linked to the magnetic catalysis phenomenon. 
Moreover, any model with the same critical behavior for a would lead to a similar effect. This means that our results 
are indeed model independent, since any relativistic theory of interacting fcrmions that belongs to the universality 
class determined by the critical behavior here considered would yield a similar kink-like feature. 

One should note that the mean field behavior of the dynamical mass a ^ \/ eB — eBc may change if higher order 
corrections (fiuctuations) are taken into account in the gap equation. The fluctuations could either change the phase 
transition to a first order one, with a discontinuity in a at the phase transition point (this was a suggestion made by 
Laughlin in |34|), or to a non-mean-field order phase transition, with the scaling law a ~ {eB — eB^Y where v > 1/2. 
While in the former case a discontinuity will appear in the thermal conductivity, in the latter case the conductivity 
will be a smooth function of the magnetic field, and a singularity will move to its higher derivatives. 



3. Low temperature limit 



Finally, we give an expression for the thermal conductivity when the temperature is much less than both F and 
eB. At low T the function cos is very sharply peaked at lu = 0, thus, expanding the rest of the 

integrand over lo and performing the integration we get 



NT r 



2+F2 

a2-F2 



(eS)2 (72 



F2 



5 

(T2-r2 

(eS)3 



3a2 - r2 
(a2+F2)3 

■2+F2 

2eB 



-F2 



eB (ct2 

}■ 



F2)2 



(42) 



It is easy to see, using the asymptotic of the ^/i— function at large values of its argument, that when B ^ the 
expression (^) goes to Eq.(|3^) in spite of the fact that we approach B = from the side B > T. For large fields we 
get 



NT 



3 ia- 



F2 



F2 



3(t2 



(a2+r2)3 eB{a^+T^y 



(43) 



Note that at large eB the thermal conductivity would approach some constant value in the case of constant mass, 
if (T > F, k{B) approaches that asymptotical value from below. This resembles the low T behavior of the thermal 
conductivity in d— wave cuprates |28| . On the other hand, if the mass a is a dynamical one with asymptotic behavior 
a ~ V eB as in our four-fermion model, then k{B) goes to zero as 1/eB at large fields. This is different from the 
thermal conductivity behavior for massless particles which tends to the universal constant k = NT/ 3 for ct = at 
large B (see, Eq.(|f)). 



III. THERMAL CONDUCTIVITY PROFILE: NUMERICAL CALCULATIONS 



In this section we do a numerical study of the profile of the thermal conductivity versus the applied magnetic field, 
taking into account the generation of the dynamical mass at a critical field that depends on the temperature. The 
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field-dependence of the finite-temperature dynamical mass is obtained from the solution of the gap equation (A4) 
derived in Appendix ^ 

To numerically investigate the behavior of the thermal conductivity ( p?] ) within a parameter range that can be of 
interest for condensed matter applications, we need to restore all the model parameters, like h, c, kB,VF,VA- Following 
Ref. pTl we write the Lagrangian density as 



(44) 



where vp and entering in v = ^vpv/^ were defined in the previous section, and the external potential is given by 
Eq.(p|). As known, this Lagrangian is equivalent to 



C 



1 



Aj) t/ij] - (tviJj^tJj^ 



Na^v 



(45) 



since the Euler-Lagrange equation for the auxiliary scalar field a obeys the constraint a — ~{g/N)^iil}i so that 
the Lagrangian density (|45| ) reproduces Eq. ( ^4|) upon application of this constraint. The effective action for the 
composite field a can be obtained by integrating over fermions in the path integral. From the minimum condition of 
the effective potential V{(j) one finds that, at fixed T, there is a critical value of the magnetic field ^/eWc/T ~ 4.1476 
such that for subcritical fields eB < eBc the gap is zero, while for eB > eBc it is given by the non-trivial solution of 
the gap equation ( |A4| ) (see Appendix ^) . 

Notice that crv has dimension of energy and plays the role of mc^ in the Dirac Lagrangian density. To generate 
a plot of k/kq (kq is the thermal conductivity at zero field) versus the magnetic field, we need to substitute T — s- 
ksT, r hT , eB — > (hv'^ / c)eB , where B is measured in Gauss. It is convenient to measure all energetic quantities 
in degrees of if, what leads to the replacement eB 0.8 • 10^^ [v / cY ° K"^ ■ B{Tesla), where the magnetic field is 
measured now in Tesla's. Thus, using the approximated value of the characteristic velocity vd — lO^cm/s ||ll[, [^ , 
eB 2.92 • 10^ • ■ B{Tesla), and we obtain the critical curve B = 0.014 • T^,. 

Let us numerically find the profile of k versus the magnetic field at a fixed T. In Fig.|l| k (B) has been plotted for two 
different temperatures taking into account the generation of the dynamical mass for eB > eBc- Here we corroborate 
what we had already argued in Section 2 based on the analytic result found for k: due to the appearance of a at 
a critical field that depends on the temperature, k exhibits a kink behavior in its profile with the magnetic field. 
Moreover, for B > Be the kink is followed by a region where k is only weakly dependent on the field (plateau-like 
region). With decreasing temperature, the position of the kink moves to the left in accordance with the critical line 
Be = 0.014T2. 
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FIG. 1. The magnetic field dependence of k at T = 20-?^ and T — 15K in the narrow width case (F — 5K). The solid lines 
represent k/kq when a QP gap a is MC-induced at B > Bc{T) (Bc(20) = 5.75T, Bc(15) = 3.23T). The dashed lines represent 
the behavior of k/kq when a remains zero at B > Bc{T). 



The numerical calculations revealed the sensitivity of the kink-plateau feature of the thermal conductivity to the 
relation between F and T. Only when F was not much smaller than T the thermal conductivity showed a kink-plateau 
profile (in Fig.|l| the curves shown correspond to the ratios F/T = 0.25 and 0.33)). 
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IV. CONCLUDING REMARKS 



In the present paper we study the thermal conductivity of relativistic fermions in a (2+l)-dimensional four-fcrmion 
interaction model as a function of the applied magnetic field, the temperature and the particle width. We have shown 
that, for certain relations among these parameters, the profile of the thermal conductivity versus the applied field 
exhibits a kink-like behavior at B ~ i?c, where Be is the critical field for the generation of a fermion dynamical mass u, 
followed by a plateau- like region at B > Be- We point out that the kink effect is the consequence of two main features: 
the generation of a fermion gap in the presence of the magnetic field (MC phenomenon) , and the enhancement of the 
zeroth-to-first LLs transitions. 

A main outcome of our investigation is that the relevant properties of the thermal conductivity of the (2+1)- 
dimensional relativistic QP system around the critical point are model independent. Indeed, the essential ingredient 
of the effective model required to produce the kink-like effect in the thermal conductivity is the critical behavior of the 
dynamical mass induced by MC near the phase transition, so not much depends on the concrete form of the effective 
Lagrangian. This fact makes our result relevant beyond the particular model under consideration, linking it to the 
universality class of theories with such a critical behavior. Such an universal character opens a window for possible 
applications. 

From a quantum field theory viewpoint, condensed matter systems whose Fermi surface are only characterized by 
nodal points are especially interesting for us, since at low energies they can be described by relativistic quantum 
field theory models of massless fermions |^^. Along that direction, a feasible possibility for the application of our 
results is the heat transport properties of graphite in the presence of a magnetic field. Let us recall that HOPG 
materials ]l9t have layered structure with two isolated points in the Brillouin zone where the dispersion is linear. 
Their electronic states can be thus described in terms of relativistic charged particles |^9| . This graphite could exhibit 
the phenomenon of MC as suggested in Ref. to explain the semimetal-insulator phase transition observed in HOPG 
in the presence of a magnetic field perpendicular to the layers. As the quasiparticles in the graphite are subjected to 
Landau quantization under a perpendicularly applied magnetic field, our results should have full strength there and 
we anticipate that the thermal conductivity of these systems will show a behavior similar to the one reported in the 
present work (for computation of the electric conductivity in graphite along the lines followed in the present work, 
see recent paper ]3T|). 

On the other hand, the characteristic feature of d-wave superconductors is also the existence of nodal points (four in 
this case) where the order parameter vanishes, thus the Fermi surface consists of four isolated points with excitations 
around them being well-defined gapless quasiparticles (QP). The kinetic part of the QPs effective Lagrangian is 
nothing but the Dirac Lagrangian for two species of massless four-component spinors | |3^ . At low temperatures such 
QPs give the main contribution to thermodynamic and transport properties. There is now considerable experimental 
evidence for the existence of well-defined QPs in the superconducting state of cuprates (see Ref. and references 
therein) . 

As was mentioned in the Introduction, the MC has been suggested ||]— |jll| to be behind the odd behavior of the 
thermal conductivity of high-Tc superconducting cuprates in a magnetic field observed in the experiment , [p^ , 
[ p^ , although alternative solutions, not based on MC, has also been proposed |24|, [Q. According to the experiment 
done by Krishana et al. [Q, and later reproduced by other authors ||l5|, ||l^, at temperatures significantly lower than 
Tc of superconductivity, the thermal conductivity k{B), as a function of a magnetic field perpendicularly applied to the 
cuprate planes of the samples, displays a sharp break in its slope (kink-like behavior) at a transition field 5^, followed 
by a plateau region in which it ceases to change with increasing field up to the highest attainable fields ~ 14T. The 
critical temperature for the appearance of the kink-like behavior scales with the magnetic field as ~ ^/B similar 
to the scaling of the critical temperature with the field found in NJL models It is worthy to mention here that 
the reliability of the thermal conductivity experiments of Krishana et al. |l4| has been a matter of debate in the last 
years ||l^] , although it seems that finally the contradictory results have been clarified and understood |l^ . 

We emphasize that none of the previously mentioned ||ri|, ||3^ attempted explanations of Krishana's 

experiment were able to obtain on a theoretical basis the kink-like behavior observed in the experiment [|l4| . In 
the case of our results, although it is striking that, as discussed in our previous paper and shown in Fig!|l|, our 
numerical curves exhibit a profile of the thermal conductivity with the applied magnetic field with similar qualitative 
characteristics!^ to those observed in the experiment , we cannot claim that our results are valid in the vortex state 



^In addition to the kink-plateau behavior and the scaling of the critical field with the temperature, our curves have another 
similarity with the experimental behavior reported in , namely, with decreasing T the crossing of the curves occurs in such 
a way that the lower T curve reaches the higher value at large fields. 
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of the superconductor since our calculations are based on the Landau level quantization, whose a ppl icability to the 
d-wave superconductor in the presence of vortices has been recently subjected to intense criticism |35-37|. 

Nevertheless, the mechanism of generating a kink-like effect in the thermal conductivity via the MC phenomenon, 
as studied in the previous sections, should be relevant for condensed matter systems with Dirac-like charged QP on 
which Landau level quantization is feasible. In this direction, we would like to point out that there is still some 
chance for the realization of MC in cuprate systems. As it has been recently argued in Refs. [Q, jso), the physical 
picture underlying the description of superconducting cuprates may involve the interplay of two different phases with 
their corresponding order parameters: one dx^_y2 superconducting (DSC) and one of density wave order (DDW). 
The DDW state, unlike the DSC state, does not break gauge invariance, so the quasiparticles of the DDW state can 
form Landau levels under an applied magnetic field ||39| . Although the above scenario is still speculative and more 
experimental confirmation is required, one can venture that given that the charged QP excitations around the nodes 
of the DDW order parameter have Landau levels, they can be subjected to a MC phase transition at some critical 
magnetic field, and hence our theoretical results could be relevant in that case for the description of the heat transport 
properties of the DDW-state in the magnetic field. 
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APPENDIX A: DYNAMICAL MASS SCALING NEAR THE PHASE TRANSITION POINT 



The effective potential of the (2-|-l)-dimensional NJL model in a constant external magnetic field at finite temper- 
ature was computed in Q . The integration over fermion fields in the functional integral with the Lagrangian (Q) can 
be performed to give the effective potential in the form 



V{a) 



NeB 



dt 



cotheSie 



64(01- 



where the temperature independent part of the potential is given by 



Vo,B{<y) + VT.B{(T), 



(Al) 



Vo.b{'^) 



Na^ NeB dt 
N 

TT 



cotheSt 



1 a' 



1 - 



aeB 



(A2) 



(the mass scale parameter Mq — h/ g — A/y^ and A is the ultraviolet cutoff which is taken much bigger than all other 
parameters in the model) and the part depending on temperature is 



Vt.b[(t) 



NeB 



dt 



47r3/2 7o t3/2 



cotheSi 



04 0| 



' A-kTH 



(A3) 



Here 04(f |t) is the Jacobi's elliptic function. 



The gap equation dV{a)/da = follows from (Al) 
Mo 



1 



/eB , cr 

tanh — 

2(7 2T 



dt e "B 



TTt e 



2t 



1 



Oi 0| 



ieB 
' AttTH J 



= 0. 



(A4) 



The solution of this equation defines the fermion dynamical mass. At T = 0, S = the gap equation admits a nontrivial 
solution only if the coupling g is supercritical, g > gc — 7r3/2/A (A/o < 0). The gap equation at T = 0, i? 7^ was also 
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studied in the literature (see, for example, Refs. It was shown that it always has a nontrivial solution at all 

g > no matter how small the magnetic field B might be. In the weak coupling phase the field-induced dynamical 
mass at zero temperature behaves as o-q = eB/2Mo in weak fields {VeB <C Mq), while at high fields {sfeB ^ Mq) 
it is (To — 0.446\/eB. At finite temperature the critical line in the B — T plane separating the massless and massive 
phases was calculated numerically in (see also We shall derive here an analytical solution of the gap equation 

near such a critical line. 

We start, actually, with the derivation of the Landau-Ginzburg-like potential by expanding V{a) in powers of cr, 
since near the phase transition point a is small, and write 

V{a) = V{0) - ^M{T,eB)a^ + ^X{T,eB)<T'^. 

As it is accustomed for second-order phase transitions, the equality of the coefficient M{T, eB) to zero defines the 
phase transition curve. The region of parameters T, B where M{T, eB) > corresponds to the spontaneously broken 
phase with fcrmions acquiring the mass, whereas the region M(T, eB) < corresponds to the massless phase. Our 
goal is to obtain the coefficients M{T,eB), and X{T,eB). When it is done the solution of the gap equation is given 
by 



M{T, eB) 
X(Tc,eB) 



(A5) 



From (A5) we obtain the behavior of the dynamical mass near the phase transition point {M{Tc, eBc) = 0). In 
particular, at fixed temperature and B — > eBc we find 



a^{T = Tc,eB) 



M{n, eB) _ _d_ f M{Tc,eB) 
X{Tc,eB) ~ deB \ X(Tc,eB) 
1 d 



{eB - eBc 



X{Tc, eBc) deB 



B=B^ 

M(Tc,eB)\cB=eB^ {eB-eBc). 



(A6) 



The potential Vo,_b(o') is easily expanded in powers of . Let us turn to VT,s(cr) which we write as the sum of two 
terms 



Vt.b{(^) 



NeB 
47r3/2 



dt 

pj2 



84(01 



-t- 



NeB 
47r3/2 



°° df 2 

-^e^'" (cotheBt-l) 

fi/2 



64(01 



1 



AttTH 



1 



(A7) 



The second term in the last expression can be expanded in a series in cr^ since the two brackets (with cotangent and 
0— function) regularize the behavior of the integrand at infinity and zero, respectively. Thus we need to calculate the 
first term 



dt 



84(01 



AttTH 



- 1 



n=l 



°° Ziyi 
4V22V^^^ifi/2 



^) =-4V^riog(i 



dt 



(A8) 



where K^{z) is a modified Bessel function {Ki^2iz) = (7r/2z)^/2g 
Finally, we obtain the following expressions for the coefficients 



n I 



M, 



Jo 



/eB 



C(l/2) 
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^ ' ' 4x/2 48 



/eB 



/eB 



'eB 



4T 7o t^' 
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{eB) 



3/2 



27r 



I dtt^/^ (cotheSt- 1) 
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(A9) 
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The critical curve M{T, B) = can be analyzed analytically at T ^ \/ eB where the integral in Eq.(A9) is exponen- 
tially small and for Tc we get the equation 



VeB 



Mo C(l/2) 
V2 ■ 



^eB 



For the solution to exist, it must satisfy that yeB ^ Mq, what gives the critical temperature 



eB 



1 

2^0, 



(All) 



(A12) 



where ctq = cr(T = 0) is the dynamical mass at zero temperature. 

One can convince oneself that there is no solution of the equation M{T, B) — when Tc ^ V eB. Indeed, for that 
let us write the integral in (RS) as 



dt 

jrj2 



(cotheSi- 1) 



64(01 



/eB 



°° dt , ^ 

^(coth<-l) 



04 0| 



4ttTH 
ieB 



- 1 



- 1 



(A13) 



We further divide the integrand into three pieces 



1 



/eB 
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coth t 
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ieB 
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Jij2 



coth t 
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-1 O4 0| 
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' AttTH 



(A14) 



The first and second integrals in the last expression can be evaluated exactly (after changing the variable i — > a;^ in 
the first integral) with the help of the formulas 



°° dt_ 



o 



^(o| 



ieB 



-4 log 2 



T 



eB 



(A15) 
(A16) 



The integral with O4— function is calculated using the Jacobi imaginary transformation to O2— function and keeping 
in it only the first term in the series when eB — )■ 0: 



O4 



ieB 



This reduces the third integral to 



eB 
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(A17) 



(when proceeding to the last equality we changed the variable t to eBt and then expanded over eB). 
Thus, combining all formulas we get the following expression for the M— function: 



M(r, B) 



NVeB 



Mq 

xfiB 



- 21og2 



/eB 



AT 



eB ^2 T 



As seen, there is no solution as eB for A/q > 0. (In case eB = and Afo < we get the standard expression for 
the critical temperature T^. = \Mq\/2 log 2 ||2^). Hence, we arrive at the conclusion that the only remaining possibility 
is that the root of the equation M{T, B) — is of the order of T^ — V eB. However, in this case we cannot expand 
the integral in Eq.(|A9|) and should turn to a numerical calculation. The root of the function M{T,B) when the 
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parameter Mq ~ is found to be \J eB/T ~ 4.1476 what defines the critical line. For the critical temperature this 
gives Tc ~ 0.54(To in agreement with the result of Ref. 1^ . We calculated numerically the coefficient before eB — eBc 
in Eq. ( |A6| ) which is found to be 0.2738, thus the scaling of the dynamical mass near the critical line is given by the 
formula 

cr ~ 0.523v/eB - eSc. (A18) 



APPENDIX B: KUBO FORMULA 



For the sake of completeness we derive here the expression for the thermal conductivity in the two-dimensional case 
used in Section II. We start from the Kubo's formula for the thermal conductivity tensor 

^^^A^) J "^^J dXTT{pUJ{0)U^it + i\)}e-''^\ (Bl) 



where V is the space volume, T the absolute temperature, p is the density matrix, and Ui are the heat current 
operators with 

U,{t) = e'"'U,e-'"K (B2) 



Integrating over t by parts in Eq.(Bl), and taking into account that the currents go to zero at t oo, we obtain 



oo " 



1 /■ e"''^* - 1 /■ d 
««.M ^yfj dt J dX-TT{pU,iO)m + m 



oo 

1 f e"*"* - 1 

^VTj '^^ ^ TT{pU,mum - US + z/3)]}, (B3) 



where we used also the fact that the quantity under Tr is a function only of i + zA. 
Now, taking into account that 

Tr{pU,{0)U^{t + «/3)} = Tr{|e-'3«[/j(0)e^^*"^^C/,(0)e-*^*+^^} = Tr{pt/,(i);7,(0)}, (B4) 



we obtain from (|B3|) and (B4) the known expression 



oo 



- -TPf / dt Tr{p[C/,(i), C/,(0)]}. (B5) 



— / dt 

VT j w 



Using that 

Tr{pU^{t)MV = Tr{C/j(0)C//(t)pt} = Tr{pC/,(0)C/.(t)} (B6) 

we have 

Tr{p[C/,(<), UjiQ)]} = Tr{pC/,(t)C/j(0)} - Tr{pC/,(i)i7j(0)}^ = 2iImTr{pC/,(t)C/, (0)}. (B7) 
Then, we can write ( p5| ) as 

oo 

2i f e^'"* — 1 

'^^Ji^) = y IniTr{pt/,(t)(7,(0)}. (B8) 



The thermal conductivity for an isotropic system is given by k = K,ii{0)/d where th e su mmation over repeated indices 
is understood {d is the number of space dimensions, in our case d~ 2). Equation (B8) then takes the form 
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(B9) 



which is equivalent to Eq.(p[), or, 
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I n >— e I n > we can write 



In the representation of the Hamiltonian eigenfunctions, e 

Tr{pC/,(i)[/,(0)} = J2 ^{e'"'+^"U,e-'"' \ n >< n \U, \ m X m\} 

n,m 

where the hermiticity of the heat current operators was used. Shnilarly 
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Now using Eqs.( Bll ) and ( B12 ) and the symmetry of the matrix elements under the interchange n m we can write 
the correlator function in the form 



g{t) = Tr{pmt),U,m} = }-Y^e-^E^+^^E^-E^)t _ g-/3(s„-iJ„)j |< „ I c/^ I ™ >|2 
The retarded Fourier transform of Q (t) is given by 
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To obtain the spectral representation of Q{fl) wc define 

e{n) = lim [g{n + «/) - g{n - ir])] 

?7— >0 
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From ( B14 ) and ( B15 ) we can write 
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dv 



V — Vl — IT] 



(B16) 



We can use this representation in order to express the thermal conductivity in terms of the spectral density 0(^2), 
which is an important step for the calculations in the Green's function formalism. Using the inverse Fourier transform, 
we can write 



oo oo oo 

g[t) = I g{n)e-'^'dn = lim J- [ e-'^'dfl [ e(iy) 
J v^o 2in J J 



dv 



h' — ~ irj' 



t > 0, 



(B17) 
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and considering the integration formula 



LO — irj 



-dcu = 



2TTie-'^\ t > 



0, 



t < 



the function Q{t) can be transformed to 

oo 

1 



g{t) = Hm ^ , 
v^o 2711 



Q{v)dv 



-iQt 



V — VL — irj 

— oo — oo 

oo 

Hm / e(f7)e~*"*-''*df^, t>0. 

77^0 J 



dfl = hm - — 

»7^0 2711 



Lu — irj 



-duj 



2VT 



(B18) 



(B19) 



To express k in terms of the spectral function we substitute with (B19) in (BIO), so we get 

oo oo oo oo oo 

j dttg{t) = J dtt Yi-m j e(f7)e-*"*-''*dr2 = Hni j Q[n)dVL j dtte-'"*-''*. (B20) 



Taking into account that the spectral function is an odd function, = — 0(— t^), we can see that only the imaginary 

oo 

part of / dte~'^*-'^H remains, what leads to 



oo 



2VT an 



or equivalently, 



K = lim ^ - e{-n)] 



(B21) 



(B22) 



From the representation ( B15 ) for the spectral function it can be shown that the thermal conductivity ( ]B22| ) can be 
expressed in terms of imaginary time Green's functions. Indeed, let us introduce the following thermal Green function 



n(r) = Tr{pe"^U^{0)e-"^U^m 



and its Fourier transform 



(B23) 



(B24) 



Inserting the complete set of energy eigenstates \m), \n) we can perform the integration over r as indicated in (B24), 
to find 



Z 



2 e-(^"'-^")^ - 1 

En - Em + ii^k 



(B25) 



If we now define an analytical function n(cj) in such a way that at dis crete points uj = icuk it coincides with Hiiujk) 
and has a branch cut along the real axis, then the spectral function (B15) is related to the discontinuity of n(w) 
across the cut 



e{n) = lim [n {n + ie)-u{n- ie)] 



(B26) 



Substituting ( [B26| ) in ( |B22D we arrive to Eq. (10) of Section II. 
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